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N00014-83-K-0258  during  the  period  of  August  1,  1983  -  January  31,  1984. 
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During  this  period  we  havcT  (1)  derived  the  backsca ttering  coefficients  for  a 
two-layer  anisotropic  random  medium;  (2)  calculated  the  emissivities  from  a 
two-layer  anisotropic  random  medium  ;  and  (3)  participated  in  the  microwave 
sea  ice  measurement  program  at  the  Cold  Regions  Research  and  Engineering 
Labora  tory  (CRREL) . 

A  two- layer  anisotropic  random  medium  model  has  been  developed  to  study 
the  active  remote  sensing  of  sea  and  lake  ice.  The  dyadic  Green's  function 
for  a  two-layer  anisotropic  medium  is  used  in  conjunction  with  the  first-order 
Born  approximation  to  calculate  the  backsca ttering  coefficients.  The  random 
permittivity  fluctuation  is  also  assumed  to  be  anisotropic  and  characterized 
by  the  correlation  functions  which  are  related  to  the  shape  of  the  fluctuation 
structure.  It  is  shown  that  strong  cross-polarization  occurs  in  the  single 
scattering  process  and  is  indispensable  in  the  interpretation  of  radar 
measurements  of  sea  ice  at  different  frequencies,  polarization  and  viewing 
angles.  The  theoretical  model  is  also  shown  to  correspond  to  the  ellipsoidal 
discrete  scatterer  model,  which  enables  us  to  determine  the  relationships 


between  the  cross-correlation  and  autocorrelation  functions.  A  manuscript  has 
been  prepared  for  submission  to  a  journal  for  publication  [Appendix]. 

The  emissivity  of  a  two-layer  anisotropic  random  medium  has  been 
calculated  using  the  dyadic  Green's  function  for  a  two- layer  anisotropic 
medium  and  the  first-order  Born  approximation.  The  emissivity  is  calculated 
by  obtaining  coherent  and  incoherent  reflectivities  and  by  making  use  of 
the  relationship  e  *  1  -  r.  The  incoherent  reflectivity  is  obtained  by 
integrating  over  the  upper  hemisphere  the  bistatic  scattering  coefficients 
obtained  under  the  Bom  approximation.  The  theoretical  results  have  been  used 
to  interpret  the  passive  microwave  remote  sensing  data  from  vegetation  canopy 
which  also  show  strong  anisotropic  dependencies.  A  manuscript  is  being 
prepared  to  document  the  theory  and  the  theoretical  calculations. 

We  have  participated  in  the  winter  microwave  remote  sensing  measurements 
at  CRREL.  Several  trips  were  made  to  the  experimental  site  at  CRREL  in  order 
to  assist  and  participate  in  the  experimental  efforts.  Our  involvement  in  the 
experimental  efforts  has  provided  us  with  valuable  insights  in  the  development 
of  theoretical  models  and  data  interpretation.  We  are  currently  waiting  for 
the  calibrated  and  reduced  data  in  order  to  start  the  interpretation  of  the 
experimental  data  with  our  theoretical  models.  / 
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Abstract 


A  two-layer  anisotropic  random  medium  model  has  been  developed  to  study  the  active 
remote  sensing  of  earth  terrain.  The  dyadic  Green's  function  for  a  two-layer  anisotropic 
medium  is  developed  and  used  in  conjunction  with  the  first  order  Born  approximation 
to  calculate  the  backscattering  coefficients.  Jt  is  shown  that  strong  cross  polarization 
occurs  in  the  single  scattering  process  and  is  indispensable  in  the  interpretation  of  radar 
measurements  of  sea  ice  at  different  frequencies,  polarization,  and  viewing  angles.  The 
effects  of  anisotropy  on  angle  responses  of  backscattreing  coefficients  are  also  illustrated. 


*  This  work  was  supported  by  the  ONR  Contract  N00014-83-K-0258,  the  NASA  Contract 
NAG5-270  and  the  NSF  Grant  ECS82-03390. 


I.  Introduction 


In  the  active  microwave  remote  sensing  of  earth  terrain,  the  random  medium  model 
has  been  used  to  account  for  the  volume  scattering  effects  of  the  terrain  media  such  as 
snow,  ice,  and  vegetation.  The  problem  of  scattering  by  a  half  space  random  medium  as 
well  as  a  layered  medium  has  been  studied  in  recent  years11'7!  .  However,  these  works 
assume  an  isotropic  constitutive  relation  whereas  the  actual  medium  may  be  anisotropic 
in  nature.  An  example  is  the  dielectric  behavior  of  sea  ice.  Due  to  the  development 
of  brine  inclusions  inside  the  ice  crystal,  it  has  been  found  that  the  dielectric  loss  of  sea 
ice  is  greater  when  the  electric  Held  is  parallel  to  the  brine  inclusions,  as  compared  to 
when  the  field  is  perpendicular  to  them  ,  implying  an  electrical  anisotropy.  It  has 
been  also  observed  that  the  c  -axis  of  the  crystal  structure  of  sea  ice  has  a  preferred 
azimuthal  orientation  .  The  radar  backscattering  coefficients  for  sea  ice  ‘ ,r,_ also 
strongly  suggests  a  theoretical  model  with  an  anisotropic  permittivity  tensor.  Besides 
these  observations,  there  are  many  other  experimental  data  which  assert  the  anisotropic 
dielectric  behavior  of  sea  ice  *IT_-’i;  .  The  cross- polarized  backscattering  coefficients  have 
been  calculated  with  the  Born  approximation  or  the  bilocal  approximation  I-3!  to 
second  order  with  the  isotropic  random  medium  model  in  layered  structures.  However, 
the  large  cross-polarization  components  as  measured  in  sea  ice  strongly  suggest  that  the 
effect  is  a  first  order  contribution. 

In  this  paper  we  study  the  problem  of  active  remote  sensing  with  a  two-layer  anisotropic 
random  medium  model.  The  dyadic  Greens  function  for  a  two-layer  anisotropic  medium 
is  first  obtained  and  approximated  in  the  far  field.  The  random  permittivity  fluctuation 
is  then  characterized  by  three  correlation  functions:  (i)  the  autocorrelation  between 
azimuthal  fluctuations  at  two  different  spatial  points,  (ii)  the  autocorrelation  between 
vertical  fluctuations  at  two  points,  and  (iii)  the  cross-correlation  between  azimuthal  and 
vertical  fluctuations  at  two  points.  With  the  information  about  the  shape  of  the  fluctuation 
structure  known,  the  third  is  related  to  the  first  two.  The  first  order  backscattering 
coefficients  are  calculated  with  the  Born  approximation.  The  results  are  examined  and 
interpreted,  emphasizing  the  effect  of  anisotropy  of  the  random  permittivity  on  the 
backscattering  coefficients  including  cross-polarization,  as  a  function  of  frequency,  incidence 
angle,  incidence  azimuthal  angle,  and  tilted  angle  of  optic  axis.  The  theoretical  model 
is  shown  to  correspond  to  the  ellipsoidal  discrete  scatterer  model,  which  also  enables  us 
to  determine  the  relationships  between  the  cross-correlation  and  autocorrelation  functions. 
Finally  the  theoretical  results  are  applied  to  the  interpretation  of  experimental  data  obtained 
from  sea  ice  measurements. 


2.  Formulation  of  the  Problem 


Consider  an  electromagnetic  plane  wave  with  linearK  polarized  lime-harmonic  field 

E«h(M)  =  (2. 

incident  upon  a  layer  of  an  anisotropic  random  medium  with  a  permittivity  tensor 


f,(r)  =<  f,(r)  >  +o/(r)  (2.2) 

where  <  r,(r)  >  is  the  mean  part  and  ?,/(?)  represents  the  randomly  fluctuating  part  [Fig. 
1],  For  die  statistically  homogeneous  medium,  <  i\(?)  >  will  be  a  constant  independent 
of  position.  f,/(?)  is  a  centered  random  function  of  position  and  its  ensemble  average, 

<  T\f{f)  >  ,  is  zero.  It  is  assumed  that  the  amplitude  of  is  small  compared  to 

<  Fi(r)  >  .  In  general,  both  <  ?,(f)  >  and  f , /(r)  are  taken  to  be  uniaxial  with  the  optic 
axis  tilted  off  the  z  -axis  by  some  angle.  For  example,  in  the  case  of  sea  ice,  the  brine 
inclusions  inside  ice  crystal  are  elongated  and  have  preferred  directions  which  are  tilted 
from  the  vertical  axis  |lsl  .  The  optic  axes  of  <  fi(r)  >  and  fi /(?)  are  rotated  by  ^  and 
4>f ,  respectively,  from  the  z  -axis  around  the  x  -axis  so  that  they  are  in  the  yz  -plane  as 
shown  in  Fig.  2. 

By  the  rotaton  of  coordinate  axes,  the  mean  permittivity  tensors  <  ?,(?')  >  and 

<  z,(f)  >  ,  before  and  after  tilting,  respectively,  are  obtained  as  follows. 
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where 

<  I  I /(?)  =  <l /,  < =  <  1/  <•<»•“*  V’/  +  ( izf  sill2  V7. 

'•„«3/(f)  =  <v>f(r)  =  (<  i,/  -  <  i/)<‘os  </-•/  sin  </'/,  (-•#) 

<33/(r)  =  <i/»m2V»/  +  ( \tf  cos2  i/>f. 

Tlie  layer  of  random  medium  has  boundaries  at  ;  =  o  and  =  =  -<i .  Hie  upper  region 
is  free  space  and  isotropic  with  permittivity  r0  and  the  lower  region  is  homogeneous  and 
isotropic  with  permittivity  <a .  All  three  regions  are  assumed  to  have  the  same  permeability 
• 

Tlie  total  electric  field  in  region  0  which  is  the. sum  of  incident  and  scattered  fields 
satisfies  the  following  homogeneous  vector  wave  equation 

V  X  V  X  7?o(r)  -  k'i  E0 (?)  =  0,  z>  0  (2.9) 

where  k2a  =  w2nt0.  Tlie  electric  field  in  region  1  satisfies 

V  X  V  X  Si(?)  -  w2fiii(f)  •  7j|(f)  =  0,  ~d<z<  0  (2.10) 

i.e. 

V  x  V  X  7?i(f)  -  t»*/iflm  •  F,(f)  =  *?(f)  •  Fi(f),  -</<*<  o  (2.11) 

where  Q(t)  =  w2n?lf(?)  is  treated  as  an  effective  source  distribution. 

We  can  express  the  formal  solutions  of  (2.9)  and  (2.11)  for  the  electric  fields  in  both 
regions  in  terms  of  dyadic  Green's  functions. 


#o(F)  =  STW  +  Jv  d3r{  Coi(r,  7 ,)  •  Q(rt)  •  El(7[) 
7?i(f)  =  F[°\r)  +  f  c/3rlf71i(r,F1).(7(f1).^(Fi) 

*  Vi 


(2.12) 

(2.13) 


where  the  integrations  extend  over  region  1  occupied  by  the  random  medium. 

The  superscript  zero  in  740,(f)  and  7 ?(,n,(r)  refers  to  the  unperturbed  solutions  in  the 
absence  of  die  random  fluctuation  part.  We  will  use  the  parenthesized  superscripts  (1), 
(2),  (3),  etc.  to  denote  higher-order  terms.  Tlie  dyadic  Green's  functions,  r;0i(F,fi)  and 
c7n(F,F,) ,  satisfy 


V  X  V  x  ('Oi(F,  fi)  --  w2n(o(j'0\(r,r])  =  0,  z  >  0 


(2.1‘) 


V  X  V  X  <V|t(r,r|)  -  |m  •  i  ( r  i )  =  /<V  -  ^  i ),  — «  <  ~  <  o  (-•>'>) 

with  the  appropriate  boundary  conditions  at  r  =  o  and  r  =  -<l  and  radiation  conditions 
at  5  =  ±oo  for  the  electric  fields.  Here  the  first  and  second  subscripts  refer  to  the  regions 
of  the  field  and  source  points,  respectively. 

We  can  obtain  the  solutions  to  (2.12)  and  (2.13)  by  iteration.  Repeating  the  substitution 
of  (2.13)  into  (2.12)  we  find  the  solution  for  the  total  electric  field  in  region  0.  in  the  form 
of  the  infinite  Neumann  series  12,11  . 

=  7?;;,,(r)+  f)7?inV)  (S-itt) 


where 
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7?o\r)  is  the  first-order  scattered  field  and  the  rest  is  the  higher-order  scattered  fields. 
Physically,  J?l0 n(f)  accounts  for  the  single  scattering,  hf\f)  accounts  for  the  double 
scattering,  and  so  on.  If  the  scattering  is  weak,  we  would  expect  that  the  series  in  (2.16) 
will  converge  rapidly.  Then  the  first  term  in  the  infinite  scries  expansion  will  provide  a 
good  approximation  to  the  scattered  field. 


F0(r)  =  hi °)(r)  +  Fo  (r)  as  Eq\t)  +  El0  °(r). 


(2.20) 


Now  we  impose  the  so-called  first-order  Bom  approximation  by  replacing  /^(f,)  in 
(2.12)  by  7i(,0,(?i) ,  as  shown  in  (2.20).  In  order  to  find  the  quantity  of  physical  interest 
we  form  the  absolute  square  of  Fa(f)  and  take  its  ensemble  average.  Then  the  mean 
intensity  of  the  electric  fild  in  region  0  will  become 
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where  we  used  Rc  <  7?jj °\r) - 7?ol)(r)"  >=  /er77n0)(F)-  <  7?o  V)’  >=  0  because  <  Q(ft)  >  = 
0 .  Therefore,  the  first-order  scattered  intensity  in  region  0  is  given  by 


<  |7?q(f)|2  >«  <  |Ty  >(r)|2  >=  j  jv  dhxdh2 
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H\  pressing  the  d\adic  Green’s  function,  the  lluctiiation  tensor,  and  the  unperturbed  electric 
lield  in  terms  of  components  in  Cartesian  coordinates,  the  scattered  field  intensity  becomes 
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3.  Dyadic  Green's  Function  and  Unperturbed  Electric  Field 


In  order  to  proceed  with  the  calculation  of  the  scattered  intensity,  we  first  have  to  find 
the  dyadic  Green's  function  (DGF)  for  a  two-lavered  anisotropic  (lilted-uniaxial)  medium. 
r70i(r,^).  with  the  observation  point  in  region  0  and  the  source  point  in  region  1.  We 
take  the  resultant  expression  for  the  DGF.  f;, „(?,/•') ,  with  the  observation  point  in  region 
1  and  the  source  point  in  region  0,  which  are  derived  in  Ref.  [25]  and  apply  the  following 
symmetrical  property  of  DGF  l-’8l 

Coi(r,  f  )  =  (,  l0(?',r)  I) 


where  the  superscript  T  denotes  the  transpose  of  the  matrix  or  dyadic.  Exchangin  ie 
arguments,  r  and  ¥  ,  in  (31b)  of  Ref.  [25]  and  taking  its  transpose,  we  obtain  t70l, 
as  follows. 
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We  notice  that  o(±/c?J  is  a  unit  vector  in  the  direction  of  an  electric  field  for  an  ordinary 
wave  and  e(fc^)  or  e(k\da)  is  one  for  an  extraordinary  wave  (upward  or  downward 
propagating).  The  coefficients  for  the  upward  or  downward  propagating  waves,  /l ’s  and 
D ’s  are  expressed  in  terms  of  the  half-space  reflection  and  transmission  coefficients  as 
follows. 
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- — (*x  +  fc2  +  *£*  -  "V.)[*2  +  *y  +  W  -  «V(<1  +  full)  ' 

«1  ~  «1*  t 

=  cos2  $k\{k\koz  -  klk\i)  +  sin2  ~  k\dx){kl  +  k0,k°lz) 

+  k\k\,{k\  -  fc2)}  +  cos  V- sin  ifky(kelx  +  kOM)(k°lz  -  k\dz)(k]  +  k\Mkot) 

cos2  Il>k%klk0,  -  k]k\i)  +  sin2  i>{k}(koz  -  k'd)(k2t  -  kotk\x) 

-  fc2A:tz(A:?  -  A:2)}  +  cos  V1  sin  ipky(k°lz  -  kOM)(k°z  +  A:'^)(A2  -  Ar„A:®,) 

cos2  ifk;{k‘ikOM  -  klk\$  +  sin2  *{*?(*„  -  k\l)[kl+ k„k'la) 

+  kyk°x(k*  —  A:")}  +  cos  i(>  sin  il>ky(k°z  +  koz)(k°z  —  k\z)(kp  +  kozk\z) 

cos2  4’k;{k\k.lx  +  kU?,)  +  sin2  u>{k'(kiz  +  Ar“)(A:2  +  k2zk°lz) 

+  klk°M{k'i  -  A:.])}  -  cos  V’sin  vky(k°z  +  kiz)(k°iz  +  A:^)(A:2  +  k°uk2z) 

:  cos2  yfk2(kjk2z  +  k2k\^)  +  sin”  v{k\(k2z  +  k\*)(kz  —  k2zk\z) 

~  klk\,(k 2  -  kl)}  +  cos  v*  Sin  tffcy(*a,  -  fc?,)(*7«  ~  ~  *a.*T.) 


r 


(3-10  j) 


lle  =  ,-os-  +  kU\i)  +  sin-  ii’{krl[biM  +  k]i)(k‘i  +  k-iakia) 

+  —  ^i)}  —  ‘‘‘>s  V’  s'[l  V’ky(k°x  +  k-MM)(k°it  +  k]dz)(k~  +■  k°tk->z) 


Under  the  far-field  approximation,  the  integral  for  (;0i(f,f')  in  (3.2)  is  evaluated  with  the 
method  of  steepest  descent  t-7-2*l  as  r  -><x  .  The  results  are 


??01(F,F')=  _ 
17 rr 


where 


Tj(k„z')  «MJfc„)LoWo(fc?,e-<^''+/>o//o(-^,)ertT.«' 


■F  Vo(kOM) 


w>l0vr «(*?,)< +  //.l^ ** 


+  .-Wc(  *;>-*-*'  +  tf.^(*£)«-<'' 
and  all  the  coefTicients  and  the  vectors  are  defined  in  (3.3)-(3.10)  with 


kT  =  k„  sin  0.  cos  <6. 


(:ui«) 


(3.116) 


/  r>  1  n  _\ 


fcy  =  k0  sin  0,  sin  <j> M  (3.126) 

and  {0 „<!>,)  are  the  observation  angles  in  the  scattered  direction. 

Given  the  wave  incident  form  region  0  with  TE  or  TM  polarization  and  the  amplitude 
Eoi ,  the  unperturbed  electric  fields  in  region  1,  7s(,0)(f)  •  are  obtained  in  a  similar  manner 
as  (7iq(f,  f/)  was. 


where 


E(t°V)TJt  = 

+  *H«i(k?'n)c'k'"'  +  BHnc(k]iyk'^y^ 

A °V)™  = 

=  =  ^t»>  ky  —  ky ,) 

Bn  ft  =  B$~(kT  —  kTil  ky  —  kyi ) 


(3.13a) 


(3.136) 


(3.11a) 
(3.1 16) 


!i  =  11  or  V 


7  =  0  or 

kp\  —  f'zli  +  kyty  (3.15a 

=  \fei~  (;}.ir,6 

{ t;i; } = ±  -  ...a,*!  -  .i.,.ty,/s  (Ms.: 

kn  =  \A«  +  i“t-  (3.I.W) 

fci.  =  sin  ^Qi  cos  0Ot  (3. 1 5c) 

*»»  =  k»  si n  0O,  si n  0Ol  (3. 15/) 

and  (0O,,  00.)  are  the  observation  angles  in  the  incident  direction. 


I 


V  *.  V 


>..vv  > 
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4.  Correlation  Function  of  Permittivity  Fluctuation 


'Hie  permittivity  nuctuation  tensor  Q(r)  =  «/*/// 1  f(r)  of  the  random  medium  is  assumed 
to  be  not  only  uniaxial  but  also  tilted  by  some  angle  wf  as  was  discussed  in  Section  2. 
The  fiiictuation  tensors  Q0(f')  and  Q(r) ,  before  and  after  tilting,  respectively,  are  given 
by 

Q(?')  0  0 

$„(?')=  0  Q(?)  0  I  (•!.!«) 

0  0  Q2(f) 

’  <?..(?)  0  0 

<?(?)=  0  QT2(f)  Qn(f)  |  ('U6) 

-  0  QT2(r)  QM 


where 


Q  u  =  Q,  Qi 2  =  Q  cos'2  0/  +  Q,  sin  2  ipf 

Qi a  =  (2:«2  =  (<2*  -  ^2)  cos  sin  V’/ 

Q:»a  =  Q  sin2  ij>f  +Q,  cos2  V/ 


(1.2) 


following  the  equations  (2.6)-(2.8). 

In  order  to  calculate  the  scattered  intensity  in  (2.23)  we  have  to  determine  the 
correlation  function  <  Qjk(7i)Q'im(f-i)  >  .  Because  the  random  medium  is  assumed  to  be 
statistically  homogeneous,  the  correlation  function  depends  only  upon  the  separation  of 
two  spatial  points  r{  and  f^291  , 


<  >=  Cjkim[7 .  -  Fa)  (4.3) 

Next  we  express  C:kim(r j  -  f2)  in  terms  of  its  own  Fourier  transform 

Cjkimifi  -  fi)  =  k'l  1 43^Wm(^)«-*?(?r,-r*)  (4.4) 

where  k\  =  Re(kx)  and  4>jfcim(/J)  is  the  power  spectral  density  of  the  correlation  function.  I 

We  define  three  correlation  functions  and  their  spectral  densities  : 

<  Q(^i)Q*(^)  >  =  Cx(F,  -  f2)  =  k\ '  j  d*n  <l>,(?)  (4.5a) 

<  Q,(7i )Q‘M  >=  ^{Fx  -  f2)  =  **  |  d*p  (4.56) 

<  Q{7 i)Q;(?2)  >  =  C3(r,  -  r2)  =  k"  J  d3/?  ♦,(3)«^'-^  (4.5c) 


where  C,(f,-?2)  and  C2(f,-f2)  arc  the  autocorrelation  functions  of  the  random  physical 
quantics  Q(?)  and  Qt[f),  respectively,  and  C2(Fi-f2)  is  the  crosscorrelation  function  of 


the  above  two  random  qualities  at  two  different  spatial  points.  r,| /-,  -  r2)  is.  in  general, 
independent  of  r,(r,  -  r->)  and  r2(r,  -  r2),  but  given  some  information  about  the  shape 

of  fluctuation  it  will  depend  on  c,  and  We  will  specify  r.t  in  terms  of  c,  and  c2 

comparing  with  the  result  of  the  ellipsoidal  discrete  scatterer  model  in  Section  6.  Now  we 
can  express  all  spectral  densities  4»jtim(,1)  in  (4.4)  in  terms  of  <l»2(,l),  and  ‘M/fy 

defined  in  (4.5)  using  (4.2). 

‘1*11 1 1  =  4*1 

‘t>i  i as*  =  4>aa 1 1  =  <t>j  <-os2  0/  +  <t>:j  sin2  ipf 

‘1*1  |;u  =  ‘1*331  I  =  ‘1*1  sin"  0/  +  ‘1*3  COS"  !p  f 

4*2222  =  4*1  ‘‘os'  0/  +  4*2  sin1  (/••/  +  24*3  cos2  ipf  sin2  V’/ 

4*22:13  —  4*3322  =  (4*i  +  4>2).siiiJ  ipf  cos2  ipf  +  <t>;t{si n 1  ipf  +  cos'*  ipf) 

‘1*3333  =  4*1  sin4  ipf  +  4><_»  cos4  pf  +  24*3  cos2  ijtf  sin2  ipf 


4*1123  =  4*1132  =  4*2311  =‘1*3211  =  (4*3  4*1 )  sill  0/  COS  Ipf 

‘1*2223  =  ‘1*2232  =  4*2322  =  4*3222 

=  [(4*3  —  4*| )  cos2  ij>f  +  (4*2  —  ‘t*:i )  sin2  rpf]  sin  ipf  cos  ipf 
4*3323  =  4*3.332  =  ‘1*2333  =  4*3233 

=  [(‘1*3  —  4>t)sin2  0/  4-  (4>o  —  4*3)  cos2  0/]  sin  ipf  cos  ipf 
4*2323  =  4*2332  =  ‘1*3223  =  4*3232  =  (4*1  +  4*2  ~  24*3)  sin2  Ipf  COS2  Ipf 

We  characterize  the  correlation  function  by  considering  two  quantities,  variance  and 
correlation  length.  Let’s  consider  the  correlation  function  of  the  fluctuation  which  is 
Gaussian  in  lateral  direction  and  exponential  in  vertical  direction. 


\p\  -  Pi? 


c,( *,  -  f,)  =  S,k"  exp  [-1*1  exp  Ltl—Sil  (1.7) 

C2{fi  -  r2)  =  Snk'f  exp  —  fy-Tj- ■  exp  -IfL—iii  (4.8) 

where  <*i  and  62  are  the  normalized  variances  of  the  lateral  and  vertical  components  of 
the  permittivity  fluctuation  tensor,  respectively,  and  l„  and  lM  are  the  correlation  lengths 
in  lateral  and  vertical  directions,  respectively. 

Here  we  took  the  same  correlation  lengths  for  two  different  correlation  functions  in 
(4.7)  and  (4.8).  It  is  physically  reasonable  to  assume  the  correlation  lengths  to  be  the  same 
no  matter  what  the  physical  quantity  is.  In  other  words,  the  correlation  length  depends 
only  on  the  randomness  of  the  medium  and  it  does  not  depend  on  the  physical  quantity 
( Q  and  Qz ).  But  it  may  differ  from  direction  to  direction  so  that  we  have  two  correlation 
lengths  in  two  different  directions,  lateral  and  vertical. 

On  the  other  hand,  the  variance  which  is  the  strength  of  fluctuation  can  depend  on 
the  physical  quantity.  <5|  ^  S2  describes  the  anisotropy  of  the  permittivity  fluctuation  of 
the  random  medium  in  our  model  and  leads  to  the  first  order  depolarized  backscattering. 
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5.  HackscaUcring  Coefficients 

Substituting  the  resultant  expressions  for  the  unperturbed  electric  field  /■;<l‘"(r)  in 
(3.13),  the  lur-field  approximated  dyadic  Green's  function  r,)  in  (3.11),  and  the 
correlation  function  defined  b>  (4.3).  into  (2.23),  we  obtain 


<«>= 


)  /  /  dz\dz->  -i  )l'k(kpt> -i  )</ti(^pi -j)^1  r-) 


where 


M°V|)  =  /  =  1,2 

I1  k[k pi i  ~()  =  ifc  •  l‘  (kpt  i  “i)>  =  1  >  2, 3 

Hi]  {kpi  -/ )  =  A  ■  £/(  ^ ’p >  )  ■  iy »  J  =  t  *  2, 3 

i,  =  x,  z2  =  if,  x3  =  x. 


(5. 2a) 
(5.26) 
(5.2.) 
(5.2d) 


Making  use  of  (4.4), 


<  |7?o(^)l2  >=  rf3? / / d2hd%e^-~k-^^  e-'« 


5 //>•* 


:2  K,(*i )**(*« )9'MF-m[z2)  *:klm((i)e 


where  p±  =  0xx  +  pvy  and  *yWm(/*)  is  given  by  (4.6).  After  the  ^-integration, 
<  |Eo(f)|2  >=  ^7jT~  /  d%  f  ^ 

El  d0*f  f 

•T*  ^ ^ 


-•/*.(»! -*a) 


The  delta  function  enables  the  /?±  -integration  to  be  performed  and  then  the  ?2  - 
integration  will  give  the  illuminated  area  A  ; 

<  pw  >-  '"trty  E  r  «.  /  / 

',r”  •/-°°  '  •'-<*  (5.5) 

l.m 


ror  ihe  case  of  backscattered  direction  where  k 
have  to  evaluate  the  integral  of  the  following  form. 


where 


C: :)- 


Jc,  —  it0 
—  *\2t 

let  —  —lc° 
K  I  xx  -  K  1  2 

*i«  =  *£,• 


The  0,  -integration  will  be  performed  term  by  term  using  the  contour  integration  of  the 
residue  calculus  and  we  pick  up  the  dominant  pole-contributions  under  some  reasonable 
approximations.  The  value  of  integration  is  significant  only  when  p  =  a,  q  =  t  or  p  = 
t,  q  —  h  and  the  result  is  the  following. 

]  ^  J  =  k'p„  +  k\\ for  p  =  .s,  7  =  f or  p  =  t,  <7  =  a  (5.8) 

l  27T«t >(0a  =  o)d,  only  for  kpin  =  (5.9) 

where  /  denotes  the  real  part  and  //  denotes  the  imaginary  part  of  that  quantity. 

Now  given  the  spectral  density  of  the  correlation  function,  we  can  determine  the  first 
order  scattered  intensity  by  a  two-layered  anisotropic  random  medium. 

The  assumptions  we  made  in  evaluating  the  above  integrals  are  given  below. 


k\m\  Kl Mi 


(5.10a) 


(The  medium  is  slightly  lossy.) 


d  » l. 


(5.106) 


(The  random  medium  layer  of  interest  contains  many  scattering  inhomogeneities  along  the 
vertical  direction.) 
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(The  scattering  layer  contains  many  wavelengths.) 


All  the  above  assumptions  are  reasonable  in  the  microwave  remote  sensing  of  the  earth 
terrain.  For  example,  in  the  case  of  sea  ice,  a  ty  pical  correlation  length  is  on  the  order 
of  millimeters  while  typical  layer  thickness  is  on  the  order  of  meters i:10i  .  Moreover,  we 
assumed  that  *t»(3)  vanishes  everywhere  on  the  hemi-circle  at  infinity.  This  assumption 
of  course  is  not  satisfied  by  all  correlation  functions.  However,  it  does  include  many 
physically  interesting  correlation  functions. 

In  addition,  because  of  die  tilted  uniaxial  fluctuation,  we  need  to  have  the  coordinate 
transformations  in  order  to  get  correct  expressions  for  «i>(/7)  which  we  use  in  (5.3)  —  (5.5). 
Since  the  fluctuation  is  uniaxial  in  x'y'z'  -coordinate  as  shown  in  Fig.  2,  the  correlation 
function  C(r,  -  ?.,)  is  expressed  as 


where 


C(ll')  =  k\A 


—in  •  It 


7 1  =  f',  —  ¥*2  =  +  yy'  +  2z' 

a  =  xaz  +  ycty  +  znt 


(5.11) 


(5.12a) 

(5.126) 


V  . 

i 

-s- 


l»y- 

& 

w 

'.V 

i 

i 

& 

r.v. 

ir.y. 

MS 


For  example,  the  correlation  function  which  is  Gaussian  laterally  and  exponential  vertically 


C(r')  =  6k 'i  exp 


p' 

2 

' 

' 

Ip 

exp 

, 

/.  . 

(5.13) 


has  the  Fourier  transform 


♦(a)  * 


61.1] 


47r2(l  +  a;/*) 


exp 


A  P 


(5.14) 


We  like  to  transform  C(^)  and  <t>(a)  into  the  form  in  xyz  -coordinate,  C[R)  and  $(p) , 
so  that  we  can  substitute  them  in  (5.3M5.5). 


C{ll')  m  C{R)  =  k\A  j  d*p4>(0)  exp  [-i/7  •  ftj 


(5.15) 


where 


7?  =  7 ,  -  =  ix  +  yy  +  zz 

if  =20  *  +  il0y  +  20. 


(5.16a) 

(5.166) 


IK 


Using  the  coordinate  transformations  [Fig.  2], 


y  =  y  cos  </>  —  z  sin  </’ 
z'  =  y  sin  0  +  z  cos  </> 

we  obtain  the  following  relations  between  a  and  fi : 

Of  z  —  fix 

<ty  —  fly  COS  0  ~  fix  SIM  0 
«,  =  fiy  sin  V’  +  fix  cos  ip 

0x  =  «z 

fiy  —  ny  cos  0  +  rt,  sin  ip 
fi2  —  —fty  sin  0  +  ft,  cos  0 


Then 


where 


dnxduydaz  =  dfixj(  y--*\lfiydfit 

\  Kz  / 

(\  I  deky  da^ 

“V'M  =  537  ^57  =  1 

/W*/  Sfc 


I  Wv  m  I 

Using  (5.17H5.20),  from  (5.11)  and  (5.15)  we  obtain 

<P(fi)  =  <t»(ax  =  fix,  c»y  =  fiy  cos  0  -  fix  sin  0,  ft,  =  fiy  sin  ip  +  fix  cos  0) 

For  the  correlation  function  given  by  (5.13),  we  have 


*(/>)  = 


exp  [-  \  filial  exp  [-\{fiy  cos  ip  -  fix  sin  0)2ffi 
Iff- [l  +  (/?y  sin  ip  +  fix  cos0)2/2] 


(5.18) 


(5.19) 


(5.20) 


(5-21) 


(5.22) 


Making  use  of  (5.8)  and  (5.9)  for  evaluation  of  the  integral  in  (5.5)  and  letting  0,  =  0O, 
and  0,  = r  7r  +  0Ol  where  0 ,  and  <p ,  are  the  observation  angles  in  the  scattered  direction, 
we  obtain  the  backscattering  coefficients  which  is  defined  by  Peake  1311  as 


OSa  —  lim 

A— « 


■l7rr2  <  l/i.'o(f)|2  >c 


fi,  ft  =  H  or  v 


(5.23) 


where  |7T0,||  is  the  incident  electric  field  intensity  with  polarization  fi,  <  |7?o(f)|2  >„ 
is  the  mean  scattered  field  intensity  with  polarization  ft  ,  and  fi  and  a  can  be  either 
horizontal  (TE)  or  vertical  (TM)  polarization.  The  final  results  are  the  following. 

1 1 

o pa  =  2 


(5.21) 
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where 


■l\  =  Moo, /'/Jo,  I* - -JTi - Iz,> -2i',°zt) 

}h  \z% 

•h  =  I AOOXH0OI  +  nao>-\(3j-<i 


(5.25/;) 


•/i  —  |  //an  //jjc 


I  -  r4** 

-i^ve 


-  i<A{k'?„,k'L,'ik'u 


t  -  c-M^.-k'r^d 

./,(  =  |/tQO,.4/3e,  +  ■4Qel4^£)l |'  - — - — Wfl(— ii"».  iu».  — i'u,  +  i'**,) 

2(*7 «  -  *i«) 
i  _  c-3(*r:.-*?::n 

•4  =  | /loo, l bn  +  ^o«,/l/3o,|“ — ^TTmI - DTT^i — ^o(-k\°2l,  ^'\zi <  _fcl«  +  *£) 

-l^iz,  "  Ki« ; 

1  _  2(k','°  +k"td)d 

Jh  =  |/WW  +  /ta«,^.,|2 — -  7~W«-.  *'.«■  +  *'.») 

-2(C,  +  *?2) 

„  i  _  e2(fcr:.+fcy;.-w 

=  1/JooA,  +  /^«/^o,|2  W,,.  *£,*£,  +  *£■) 

-lK  iz»  +  Kt*i  / 


^  8  — 


I  -  *•<*', 'IN 

I2  — — —  «c(*K.  2*K!) 


-<*K? 


i  _  cW;:? +*"::)* 

•/.  =  -  .  IlcWLWr'k'li  +  ® 

“2(*1«  +  *E?) 


I  _  -'•*v:r<‘ 

^10  =  |/WWI2  -ltc(k\~,k\<:x,2k\”) 

J 1 l  =  A  ■  2/?c(/luet/?tf0,  H veiA h ol)d  R'b 

Ra[X,Y,0x)  =  ^Y^y^{C\\A$\\i\[(i,)  +  2C\xACviA(b\\ri{fl*)  +  ClzA’b'nwiP*) 

+  2CuA£33A*frii33(/?z)  +  2 (7  2aC33a*I*2233(/?z)  +  ^La^sss^z) 

+  ZC’ttAC’sSA'&im^*)  +  2C22A^23A<1,2223(^*) 

+  "CZZACiSA’blMiiPM)  +  C23A^?a23(/?*)} 

Rb(X,  Y,0m)  =  {C'*ls4>nu(/3,)  +  2C^\ bCt’b^wtA^*)  + 

+  2CiiflC3afl<l>u33(/3«)  +  2C22XjCa3s‘{>2‘J33(/?*)  +  C'33B‘I>3333(/^*) 

+  2Ciiflf723fl‘l>u23(/3,)  +  ‘2  C22  flC'JS  O  *1*2223  (/^*) 

+  2 (^33  fl(?23B  *1*2333  (/^s)  +  C’23B<*>2323(^»)} 

Rc{X i  Y t  ftz)  —  {^UC^I  1 1 1  (P*)  +  2C|  \cC>‘iC<b \\Vl[0z)  +  ^22C<l*2222(/3*) 

+  2C’|  10^330*1*1 133(/?*)  +  2  ('•>  2  C  ^3  3  C  ‘I*  2  2  3  3  ( /^2 )  +  ('hc^ZZ-iz{0z) 

+  2(7||cC’23C‘I>II23(/^z)  +  2  Ci<c('‘Z:\  C  *1*2223  (At  ) 

+  2(733cC,23C*I*2333(/^z)  +  ^33C*I>,<M23(A«)} 


(5.25»i) 


(5.25c) 


(5.25/) 


(5.25ff) 


(5.25A) 


(5.25t) 


(5.25/) 


(5.25k) 


(5.26a) 


(5.266) 


(5.26c) 


if'  _ _ : 

"  “  1 

T.t{—  1  l  (ti)  +  ~  “‘I’l  l-J^(O)]  +  (A-xifcpt)*'l>j;{-.>:i(0)} 

ir 

(5.26ii) 

ff(A')  =  {A2,  4-  (Ay,  cos  V  -  v  sill  (/')'} 1/2 

(5.27a) 

r  , 

( '(x )  =  — '—(k;t  +  .V  •  -  «v  .)[*;,  +  A'  *  —  ,JtL{< ,  + ,  „)j 
i f i  f i*  j 

(5.276) 

r1M(.V,K)  =  (Ayjcos0  —  X  sin0)(Ay,  cos  0  -  K  sin0) 

(5.28a) 

r.joyt  =  ^  <:os‘!  ^ 

(5.286) 

=  A2,  sin2  0 

(5.28c) 

^■'23A  =  -2A2,  sin  0  cos  0 

(5.28(f) 

f’l  K)  =  Axt(Ay,  cos 0  -  A'  sin  0)(Ay,  sin  </’  4-  K  cos  0) 

(5.20a) 

C-i2o(X ,  Y)  =  — AIt  cos  V»[Ay»( Aryi  sin  V’  +  V  cos  c)  -  Aj  sin  ip\ 

(5.206) 

('33b(X,  Y)  =  A*t  sin  0[V(Ay,sin0  4-  V'  cos0)  -  Af  cosi/>] 

(5.20c) 

(•23b{X i  K)  =  Az,(A2t  sin2  ip  -  Y 2  cos2  ip  +  Af  cos  2^] 

(5.20(f) 

Ciic(A\  V)  =  A2,(Ayj  sin  0  4-  A'  cos0)(Ayi  sin  0  +  K  cosv) 

(5.30a) 

c»c(*,n 

=  (Ayi{Ayi  sin  0  4-  AT  cos  0)  —  A2  sin  0][Ay,(Ay,  sin  0  +  V  cos  0)  -  A2  sin  «-•) 

(5.306) 

C33C(X,Y) 

=  [X(Ay,  sin  0  4"  X  cos0)  —  A2  cos  0][)'(Ay,  sin  ip  4-  Y  cos  0)  —  k\  cos  ip] 

(5.30c) 

Cnc(X,Y) 

—  [Ay, (Ayt  sin  0  4-  X  cos  0)  —  Aj  sin  0j[^(Ay,  sin  0  4-  Y  cos  ip)  -  A2  cos  v) 

4-  [X  (Av,  sin  0  4-  X  cos  0)  —  A2  cos0][Ay,(Ay,  sin  ip  4-  Y  cos0)  -  A2  sin  rj 

(5.30<f) 

4>j  *im(/?*)  =  <b]klm(73_ l  =  -^>1,0*) 

(5  31 

A 

( l,  only  for  apa  =  ohv  or  ov h  when  0  =  0,  <  1  =  <  1, 

1 0,  otherwise 
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6.  Discussion  ami  Application  of  (he  Results 

Wo  have  obtained  the  backscattering  coellicients  for  a  two-layer  anisotropic  random 
medium  in  (5.24)-(5.32).  In  the  limiting  case  when  both  mean  and  fluctuation  are  isotropic, 
i.e.  when  <u  =  <u  (uf  =  ,  lf,  and  =  <i’f  =  o  ,  the  rtjults  reduce  to 


+  Hd  |  It ,  2 » | ' '  e  " 1  k • d' $  j2 1 1 

4  A  1 1  + iw.-H  i*;,  + 


(6.1a) 


>'oi. 

4  k'o  J 

[  l  -  c~ik"‘>d 

I'-ii 

N8  I 

{  2^ 

+  H</|5|2,|2c  Ak”‘'d\k\n  -  ^.sin2t?0l|'‘l>(,J)| 


(6. 16) 


OtfV  =  °VH  —  0 


These  are  exactly  the  same  as  the  results  obtained  for  die  case  of  two-layer  isotropic 
random  medium  by  Zuniga  and  Kong  lfii  .  As  an  another  limit  we  let  <p  =  ()  for  the  case 
of  vertically  uniaxial  mean  and  tilted  uniaxial  fluctuation  (v>/  #  o) .  The  results  check 
with  Lee’s  [3a|  . 

Comparing  with  the  result  of  the  active  microwave  remote  sensing  of  two-layer 
homogeneous  medium  containing  ellipsoidal  discrete  scatterers  1331  for  the  case  of  prolate 
spheroid  which  seems  to  be  similar  to  our  anisotropic  continuous  random  medium  model, 
our  result  with  v  =  0  has  one-to-one  correspondence  with  them  term  by  term  except  the 
one  which  comes  from  the  coherent  effect  between  oppositely  scattered  waves.  From  the 
comparison  of  two  results  term  by  term,  we  get  the  constraint  on  the  spectral  density  of 
<  Q{r iW.fa)  >  (4.5c): 

+1(3)  =  <t>,(3)‘M3).  (6.2) 

Assuming  that  the  three  correlation  functions  in  (4.5)  have  the  same  statistical  behavior, 
(6.2)  leads  to 

*3  =  **  (6-3) 

Therefore,  given  that  the  shape  of  fluctuation  is  a  spheroidal  ellipsoid,  the  cross  correlation 
between  Q  and  Q-:  will  be  determined  by  <l>2  =  <l>t<l>2 .  We  shall  use  diis  constraint  in 
assigning  statistical  parameters  in  order  to  get  various  responses  of  a  in  the  following. 
We  shall  choose  the  correlation  function  to  be  exponential  both  laterally  and  vertically  as 

follows. 

~  _ f  1*1  \y\-Vi\  1*1 -Sail 


CX(T\  -  f-i)  =  6\  k>Aexm  — 


(6.1a) 


The  spectral  densities  are  given  as 


i  +  1  +  «;/;!)( I  +  n'Hi) 


J  =  1,2. 


Tliese  satisfy  the  assumptions  made  in  the  calculation  of  <x . 

Now  we  look  at  various  effects  of  anisotropy  of  random  permittivity.  First  we 
see  <jhv  =  °vh  which  agrees  with  the  Lorentz  reciprocity  theorem  since  the  medium 
is  reciprocal  and  we  consider  the  scattering  in  the  backscattering  direction.  The  most 
significant  effect  is  that  there  is  a  depolarization  in  the  first-order  backscattering.  As  we 
see  in  (6.1c),  the  first  order  single  scattering  with  the  isotropic  random  medium  model 
does  not  induce  any  depolarization  effect  in  the  backscauered  direction.  Even  for  the 
case  of  uniaxial  mean  and  fluctuation  with  their  optic  axes  in  the  c  -direction,  i.e.,  when 
(|„  f)/  ^  f,4/  but  ip  =  4>f  =  ()  (untilted),  (5.24)  again  gives  <r,{V  =  aVH  =  0  with 
some  algebraic  manipulations.  This  shows  no  depolarization.  It  is  physically  clear  that 
if  the  medium  is  vertically  (untilted)  uniaxial  either  in  mean  or  in  fluctuation  the  single 
scattering  in  the  backscattered  direction  does  not  change  the  polarization  of  the  incident 
wave.  In  other  words,  there  is  no  coupling  between  two  different  linear  polarizations. 
However,  if  either  the  mean  (background)  or  the  fluctuation  becomes  tilted  uniaxial,  i.e., 
4>  ^  o  or  4>f  o  ,  then  it  gives  us  nonzero  depolarization.  We  can  consider  three  cases 
which  give  the  cross  polarization  effect  in  the  first-order  backscattering. 

First,  ip  ^  0  and  ipf  =  o  where  the  background  permittivity  is  "anisotropic"  (tilted 
uniaxial)  and  the  permittivity  fluctuation  is  isotropic  or  untilted  uniaxial.  Then  aHV  = 
aVH  o  .  This  is  because  the  wave  incident  with  either  TE  or  TM  polarization  upon  an 
anisotropic  background  medium  produces  both  the  ordindary  wave  and  the  extraordinary 
wave.  This  is  the  so-called  double  refraction. 

Second,  ipf  ^  0  and  ip  =  0  where  the  fluctuation  is  anisotropic  and  the  background 
is  isotropic  or  untilted  uniaxial.  In  this  case  the  wave  with  one  polarization  can  produce 
the  other  polarizadon,  when  it  is  scattered  in  an  anisotropic  fluctuation,  even  in  the 
backscattered  direction. 

Third,  ip  =  ip f  ^  o  .  It  is  obvious  that  when  both  the  background  and  the  fluctuation 
are  anisotropic  there  would  be  a  strong  depolarization  effect.  It  is  physically  more 
reasonable  to  have  ip  =  4>f  than  ip  ^  ip  f , 

It  is  also  to  be  noted  that  even  if  one  of  three  conditions  is  satisfied,  when  <p0i  =  ±90° , 
i.e.,  when  the  propagation  vector  and  the  optical  axis  of  mean  or  fluctuation  are  in  the  same 
plane  ( yz  -  plane),  aHV  =  () .  This  is  because  for  this  case  the  electric  field  is  polairzed 
either  perpendicularlly  (TE)  or  parallel  (TM)  with  respect  to  the  plane  of  incidence  so 


that  the  double  refraction  does  not  occur  and  the  single  scattering  does  not  change  the 
polarization  of  the  incident  wave  in  the  backscattered  direction. 

As  an  illustration  of  our  theoretical  results,  the  typical  responses  of  am,,  avv  and 
c rHV  are  shown  as  a  function  of  frequency  at  0Ql  =  ;jo°  [Fig.  3]  and  as  a  function  of 
incidence  angle  at  /  =  10  GHz  [Fig.  4].  by  choosing  the  appropriate  parameters  for 
6\,  6-2,  lp  and  /,  which  are  reasonable  numbers  in  the  micowave  remote  sensing  of  the 
earth  terrain.  We  clearly  see  that  we  obtain  a  strong  a,{V  . 

For  the  backscattering  coefficient  for  the  isotropic  half-space  random  medium  (d  —  oo) 
as  a  function  of  incidence  angle  (0o,),  avv~  is  always  greater  than  a,II{  ,  because  the 
vertically  polarized  TM  waves  are  transmitted  and  backscattered  more  than  horizontally 
polarized  TE  waves.  In  contrast  with  the  isotropic  case,  the  anisotropic  half-space  random 
medium  gives  us  the  possibility  of  having  aHH  >  avv  as  it  will  be  shown  in  Figs.  5-10. 

First,  it  is  due  to  an  anistropic  fluctuation.  In  Fig.  5,  when  >  62  ,  for  large 
incidence  azimuthal  angle  0O«  =  so°  (near  <)()° ),  aHH  is  seen  to  b'e  greater  than  avv 
especially  for  small  incidence  angles  0,„  .  This  is  because  the  strength  of  fluctuation  is 
stronger  in  lateral  directin  (plane  perpendicular  to  the  optical  axis)  compared  to  that  in 
vertical  direction  (along  the  optic  axis  of  fluctuation),  in  other  words,  the  fluctuation  is 
more  correlated  in  lateral  direction.  Near  <z>0«  =  90° ,  the  polarizatin  of  TE  wave  will  be 
almost  in  lateral  direction  with  respect  to  the  optical  axis  of  Q  .  However,  for  small  </>0i , 
it  follows  the  common  response,  i.e.,  aw  is  always  greater  than  a,n,  over  all  incidence 
angles  0Oi  as  shown  in  Fig.  4. 

In  Fig.  6,  when  <  6* ,  for  small  <t>0i  (near  0°),  a,{H  can  be  larger  than  avv 
at  small  incidence  angles  so  that  there  is  a  cross-over  between  aHH  and  avv  near 
00,  =  30° .  The  reason  for  this  effect  is  similar  as  before.  TE  waves  have  more  components 
parallel  to  the  optical  axis  near  <t>0 ,  =  0°  than  at  other  angles  and  the  fluctuation  is  more 
correlated  in  the  direction  of  optic  axis.  In  contrast,  for  large  <£0.,  <*w  is  always  larger 
than  aHH  for  all  0Oi  as  shown  in  Fig.  7.  We  note  that  this  kind  of  response  can  happen 
when  w  #  o  ,  i.e.  for  the  tilted  uniaxial  fluctuation. 

Second,  this  phenomena  also  arises  due  to  an  anisotropic  background  [i>  ^  o,  e(  ^ 
eu) .  For  example,  in  Fig.  8,  for  large  auH  is  greater  than  avv  over  a  wide  range 
of  lower  0Oi  when  t'[,  deviates  largerly  from  c"  .  See  the  case  of  =  o.oi,  t'[M  =  o.i 
compared  to  the  case  of  e"  =  o.oi,  t'{,  =  0.02 .  This  is  because  TM  wave  experiences 
greater  absorption  than  TE  waves  due  to  higher  loss  along  the  optic  axis  of  the  background 
medium. 

The  possibility  of  aHH  >  aVv  can  be  also  seen  from  the  responses  of  a  as  a  function 
of  incidence  azimuthal  angle  <j>Qx  at  a  fixed  0Ot  in  Figs.  9-10.  When  <5,  >  [Fig.  9], 
<jhh  increases  and  aYV  decreases  as  <p0l  ,  so  that  aliH  >  avv  at  large  <j>nt .  When 
<  <5a  [Fig.  10],  oHh  >  ovv  at  small  <f>0x  .  As  a  whole,  the  cross-polarization  (at{V)  is 


very  sensitive  to  0,h  ,  while  the  like-polarizations  *vv)  are  not  much  affected  by 

changing  <£«, . 

Given  the  incidence  angles,  0{) ,  and  e»„,  .  all  a 's  are  affected  greatly  by  changing 
the  tilted  angle  of  optic  axis  as  shown  in  Figs.  11-12.  We  observe  a  lot  stronger  cross¬ 
polarization  when  v  #  o  [Fig.  11]  compared  to  when  </’ =  o  [Fig.  12].  For  some  range 
of  tl’,<rVH  is  as  large  as  aHn  or  aVv  .  It  could  be  even  larger  in  some  cases.  From  this 
we  observe  that  the  anisotropy  of  background  medium  produces  more  cross-polarization 
than  the  anisotropic  fluctuation  does. 

As  an  application  we  try  to  explain  the  experimental  data  by  matching  our  theoretical 
results  with  experimental  data  collected  from  the  Arctic  sea  ice  by  Onstatt,  Moore  and 
Weeks  [l5l . 

In  Figs.  13-15,  we  have  matched  the  backscattering  data  of  the  thick  first-year  (TFY) 
sea  ice  as  a  function  of  an  incidence  angle  at  three  different  frequencies:  /  —  9  GHz 
[Fig.  13],  /  =  13  GHz  [Fig.  14],  and  /  =  17  GHz  [Fig.  15],  simultaneously.  The 
letters  //,  v  ,  and  ('■  represent  experimental  data  for  aHI{,  avv  and  aHV  ,  respectively, 
and  the  continuous  curves  represent  the  theoretical  results.  In  order  to  match  these  data 
we  assigned  a  set  of  parameters:  l„  =  imm,  lz  =  3mm,  A,  =  o.t.  62  —  0.8  with  <p&  = 
75 °,v»  =  il>f  =  25° .  The  ground  truth  for  <t> 0„  ip  and  ipf  is  not  known.  As  seen  from 
the  figures,  the  theoretical  results  fit  with  experimental  data  pretty  well  except  at  lower 
angles.  Higher  values  for  a  of  experimental  data  at  lower  incidence  angles  compared  with 
the  theoretical  values  are  due  to  the  rough  surface  elTect  of  sea  ice  surface.  Our  theory 
only  takes  into  account  the  volume  scattering  and  the  rough  surface  scattering  effect  is  not 
included.  Notice  that  we  have  chosen  /,  >  lp  and  6->  >  A,  .  Since  the  brine  inclusions 
inside  an  ice  crystal  look  like  vertically  elongated  ellipsoids  t,9l  ,  the  vertical  correlation 
length  should  be  greater  than  the  lateral  correlation  length.  The  strength  of  fluctuation  in 
vertical  direction  is  stronger  than  that  in  lateral  direction,  because  this  quantity  determines 
the  magnitude  of  scattering  due  to  the  brine  inclusions.  We  have  chosen  the  different 
permittivities  for  sea  ice  (ti,c,,)  and  sea  water  (c2)  at  different  frequencies,  following 
Vant  et  al.  t,9l  The  real  part  of  ft  is  chosen  to  be  isotropic  while  the  imaginary  part  is 
chosen  to  be  anisotropic,  following  Sackinger  and  Byrd’s  report 181 . 

In  Figs.  16-18,  we  have  also  matched  the  backscattering  data  of  the  multiyear  (MY) 
sea  ice  as  a  function  of  an  incidence  angle  at  three  different  frequencies:  /  =  9  GHz 
[Fig.  16],  /  =  13  GHz  [Fig.  17],  and  /  =  17  GHz  [Fig.  18],  simultaneously.  Another 
set  of  parameters  are  assigned  to  match  data:  lp  =  3 mm,  /,  =  t mm,  dt  =  o.i,  <52  =  0.3 
with  «t> o,  =  75°,  ip  =  4>f  =  35° .  The  multiyear  ice  has  properties  different  from  those  of 
the  first-year  ice  in  that  it  contains  mainly  air  bubbles  and  it  has  much  lower  salinity  i'91  . 
Hence,  lp  and  tz  do  not  differ  much.  The  dielectric  constant  and  the  dielectric  loss  are 
expected  to  be  less  than  those  of  the  thick  first-year  ice. 


7.  Conclusions 


Introducing  the  dyadic  Green's  function  formalism  and  applying  the  lirst  order  Born 
approixmation,  we  have  solved  the  problem  of  electromagnetic  wave  scattering  by  a  two- 
layer  anisotropic  random  medium  for  an  application  to  microwave  active  remote  sensing  of 
earth  terrain.  With  an  advantage  of  wave  approach,  we  obtained  the  analytic  expressions 
for  backscattering  coefficients  which  includes  depolarization  effect 

There  was  a  significant  result  in  that  the  anisotropy  of  mean  permittivity  or  permittivity 
fluctuation  of  random  medium  can  lead  to  nonzero  depolarization  in  the  backscattered 
direction  for  the  first  order  scattering,  in  contrast  with  the  isotropic  random  medium  which 
does  not  have  depolarization  in  the  first  order  backscattering.  Either  the  mean  permittivity 
or  the  permittivity  fluctuation  of  the  random  medium  has  to  be  at  least  uniaxial  and 
tilted  with  respect  to  vertical  (;)  direction  in  order  to  obtain  the  first  order  depolarized 
backscattering.  The  result  also  shows  the  possibility  of  < tHI{  >  avv  for  an  anisotropic 
random  medium  even  in  the  case  of  half-space.  This  does  not  occur  for  an  isotropic  half¬ 
space  random  medium.  As  an  application  to  remote  sensing,  we  have  matched  theoretical 
results  very’  well  with  experimental  data  collected  from  sea  ice. 

Since  most  of  the  experimental  data  for  backscattering  coefficients  contains  only  the 
incidence  angle  (0ot)  response,  but  not  the  incidence  azimuthal  angle  (0()t)  response,  it 
is  difficult  to  take  advantage  of  our  anisotropic  random  medium  model  at  this  point.  It 
would  be  recommended  that  some  experimental  data  for  an  azimuthal  angle  response  be 
collected  in  order  to  account  for  the  effect  of  anisotropy  of  random  permittivity. 
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Figure  Captions 


Figure  1  Scattering  geometry  of  a  two-layer  anisotropic  random  medium.. 

Figure  2  Geometrical  configuration  of  permittivity  tensor.. 

Figure  3  Theoretical  results  of  <rHH,  aVv  and  aHV  versus  incidence  angle 
(0n  =  30°,  =  op  .. 

Figure  4  Theoretical  results  of  aliH  and  «rVv  versus  incidence  angle  (0„  at  / 
=  10  GHz,  fa  =  0° .. 

Figure  5  Theoretical  results  of  oHh  and  avv  versus  inciddence  angle  with 
<t>oi  =  8tf  when  6i  >  62  .. 

Figure  6  Theoretical  results  of  oHh  and  aVv  versus  incidence  angle  with  = 
10°  when  <  62 .. 

Figure  7  Theoretical  results  of  aHH  and  <rvv  versus  incidence  angle  with  <£*  = 
l(f  when  <  Sj .. 

Figure  8  Effect  of  anisotropic  background  on  incidence  angle  response.. 

Figure  9  Theoretical  results  of  <jhh,  °w  and  aHv  versus  incidence  azimuthal 
angle  (<£*)  at  0„  =  30°  when  £1  >  S2 .. 

Figure  10  Theoretical  results  of  aHH,  avv  and  <tHv  versus  incidence  azimuthal 
angle  at  0„  =  iff5  when  <5,  <  6 2 .. 

Figure  11  Theoretical  results  of  i?hh,  °w  and  aHV  versus  tilted  angle  [p  =  rpf) 
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at  o„  =  30°  f  4>„  =  10° .. 

Figure  12  Theoretical  results  of  aHH,  <rvv  and  aHV  versus  tilted  tingle  (ipf) 
with  v  =  0 0  at  0„  =  30°,  <j>ai  =  iff* .. 

Figure  13  Interpretation  of  backscattering  data  of  thick  first-year  sea  ice  as  a 
function  of  incidence  angle  at  /  =  9  GHz.. 

Figure  14  Interpretation  of  backscattering  data  of  the  same  sea  ice  at  /  =13 
GHz.. 

Figure  15  Interpretation  of  backscattering  data  of  the  same  sea  ice  at  /  =17 
GHz.. 

Figure  16  Interpretation  of  backscattering  data  of  multiyear  sea  ice  as  a  function 
of  incidence  angle  at  /  =  9  GHz.. 

Figure  17  Interpretation  of  backscattering  data  of  the  same  sea  ice  at  /  =13 
GHz.. 

Figure  18  Interpretation  of  backscattering  data  of  the  same  sea  ice  at  /  =17 
GHz.. 
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